Spin-charge separation at small lengthscales in the 2D t — J model 
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We consider projected wavefunctions for the 2D t — J model. For various wavefunctions, including 
correlated Fermi-liquid and Luttinger-type wavefunctions we present the static charge-charge and 
spin-spin structure factors. Comparison with recent results from a high-temperature expansion by 
Putikka et al. indicates spin-charge separation at small lengthscales. 
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T—l . INTRODUCTION 

The properties of correlated electrons in two-dimensions (2D) are controversially discussed. Anderson Q has 
I/"") ■ proposed that the well established Luttinger-liquid properties of one dimensional electron systems [^j might occur, 
to a certain extend, also in 2D. In ID the spin and the charge excitations move independently at small frequencies 
[0], a phenomena called spin-charge separation. Separation of the spin and the charge degrees of freedom has some 
' interesting consequences on the properties of the ground-state wavefunction. In the framework of the t — J model 
' the charge carriers are given by the holes which form in ID, since they move independently with respect to the spin 
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rigorous in the U/t — > oo limit of the ID Hubbard model 
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background, a hole Fermi- surface distinct from the Fermi-surface of the spins. This notion has been shown to be 



A huge amount of studies have been devoted, in the last years, to the enlargement of our knowledge upon 2D 
correlated electron systems. It is probably fair to say, that at the moment there is no method in sight capable of 
I i establishing rigorously the properties of 2D correlated electron systems at large lengthscales, or, what is equivalent, 
at small differential wavevectors as it would be necessary for the determination of true spin-charge separation in 2D. 



It is, on the other hand, an interesting question to ask, whether at small to intermediate lengthscales spin and charge 
excitations can move independently or not. 

Considering a variational wavefunction originally introduced by Hellberg and Mele in the context of the ID t — J 
model H , it has been shown Q that (i) it is possible to define Luttinger liquids in two dimensions variationally and (ii) 
that the projected kinetic energy favours Luttinger-liquid type correlations in 2D, as it does in ID ||. Chen and Lee 
have shown that the static spin-spin structure factor, Sq, of this Luttinger-liquid wavefunction compares favourably 
with the Sq of the numerically obtained exact ground state of a 10 x 10 cluster. Also, the recent high-temperature 
expansion study by Putikka et al. || has provided reliable data on Sq and the static charge-charge structure factor, 
Nq, of the 2D t — J model. An unexpected result has been the observation of a prominent enhancement of Nq near 
(it, 7t) at quarter filling. Putikka et al. found additionally that Nq could be fitted very closely at all densities with the 
curve appropiate for spinless Fermions, i.e. assuming an independent hole Fermi-surface. These results for Nq have 
been confirmed at small densities by a study using the power method, Quantum Monte Carlo and a perturbation 
expansion M. 

The high-temperature expansion || and the numerical studies ||] can, of course, not determine the properties of 
the 2D t — J model at long lengthscales and therefore address the question of true spin-charge separation. They do, 
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on the other hand, provide reliable data on Sq and iV q at short to intermediate lengthscales. Here we will interpret 
these data for the correlation functions of the 2D t — J model within the context of projected wavefunctions. For a 
class of variational wavefunctions with both short- and long-ranged correlations we calculate Sq and Nq and find that 
short-ranged charge-charge repulsion can explain the data by Putikka et al. and by Chen et al. [0]^). Furthermore 
we show that it is the projected kinetic energy which gives rise to the prominent peak in A q . Interpreting these 
results we come to the conclusion that the strength of the short-ranged correlations can only be explained when the 
the spin and the charge degrees of freedom move independently at small lengthscales. These results may partially 
support Anderson's concept jl(| of incoherent one-particle tunneling in between the different CuO-layers of the high- 
temperature superconductors. 



WAVEFUNCTION 



We consider variational wavefunctions for the 2D t — J model, which is given by 

(«,i)><7 

where the cj a {t i a ) are the creation (anihilation) operators on site i of electrons with spin a =f , j, in the subspace 
of no double occupancy. The Si are the spin operators on site i and denotes pairs of n.n. on the square lattice. 
The Jastrow-Luttinger wavefunction is given by || 



|*(T,S)> = exp 



i<3 



Po |V>0>, (2) 



where |-0o) denotes the filled Fermi sea of electrons, Po the projection operator on the subspace of no double occupancy. 
hi = cj . + cj | is the density operator on site i and 

/(r,-r 3 -) = -(l-5)(l-% J - > ) In \ Ti - rj \ + S6 {itj) . (3) 



In the limit T — > oo the wavefunction defined by Eq. (|2|) reduces to the Gutzwiller wavefunction, Po \ipo). A finite 
1/(2T) > controls then the strength of the repulsive JTTJ density-density Jastrow factor, given by Eq. (||). The second 
variational parameter, S, controls the details of the Jastrow factor. For < S < 1 the density-density correlation is 
long-ranged and it has been shown ||^| that this logarithmic correlator leads to an algebraic singularity in the 
momentum distribution function at the Fermi-edge and therefore to a Luttingcr- liquid state. For S — 1, on the other 
side, the Jastrow correlator is short-ranged, with only nearest neighbor repulsion (for T > 0), and the resulting state 
is a correlated Fermi- liquid state |J . 

The Jastrow-Gutzwiller wavefunction, as defined by Eq. cannot be evaluated analytically. One then considers 
large but finite clusters and calculates properties of |\P(T, S)) numerically by the variational Monte-Carlo method 
p3| . Here we are interested in the static spin-spin structure factor, 

Sq = 1 /L e * q(r ' " Fj ' < £f 5? > , (4) 
and the static charge-charge structure factor, 
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iV q = l/L^y q ( r '- r ^ < AfiiAhj >, (5) 

where L is the number of sites of the lattice considered and Ahi = hi — n, with n being the particle density. In order to 
minimize finite-size effects in cluster calculations it is advantageous to consider closed shell configurations []l3|| . Here 
we consider clusters with antiperiodic boundary conditions which tile the square lattice in Li = (2m, 2m — 2) and 
L2 = (—2m + 2, 2m) directions (with m an integer). In Fig. (|l|) we illustrate the first Brillouin zone for a lattice with 
m = 6, which has L = 12 2 + 10 2 = 244 sites. Also shown in Fig. (|l]) are the electron Fermi-sea (filled circles) and the 
hole Fermi-sea (filled squares) for Nh = 124 holes. This type of lattices has the property that, whenever the number 
of holes is a multiple of four, but not of eight, both the electron Fermi-sea and the hole Fermi-sea are closed-shell 
configurations, as illustrated in Fig. ([[]). The dashed lines in Fig. ([j]) denote the locations of the Fermi surfaces for 
the same nominal densities n = (L — Nh)/ L of electrons and of Nh/ ' L holes respectively, in the thermodynamic limit. 
Here the Fermi surfaces in the thermodynamic limit are consistent with the respective Fermi surfaces of the finite 
cluster (this is not necessarily always the case H). 



RESULTS 



In Fig. we present for various q in the first Brillouin zone the static charge-charge structure factor, iV q , as 
defined by Eq. (^), for correlated Fermi-liquid wavefunctions, \^{T,S = 1)) as defined by Eq. (|^), for the lattice 
illustrated in Fig. (|l|) with L = 244 sites and Nh — 124 corresponding to a particle density n = (L — Nh)/L w 0.49. 
In Fig. (||) we have included for comparison the curve for iV q for spinless Fermions, i.e. assuming a well defined 
hole Fermi-sea as illustrated in Fig. ([!]). We observe that data for the Gutzwiller state, realized for 1/(2T) = 0, is 
qualitatively different from the curve for spinless Fermions. The data for |^/(1/(2T) = 0.2,5 = 1)) is, on the other- 
hand very close to the curve for spinless Fermions and agree remarkable well with the results of Putikka et al. and 
Chen et al. &@]. Noting that the correlated Fermi-liquid wavefunction |^(1/(2T) = 0.2,5 = 1)} contains explicitly 
only nearest neighbour correlations we conclude that the static density-density correlation factor of the 2D t — J model 
is dominated by short-distance effects. This does, of course, not exclude the possibility of long-range correlations to 
be present in the paramagnetic state of the 2D t — J model, as these correlations cannot be calculated rigorously with 
present-day methods. 

The energetics behind the results for iV q is dominated by the projected kinetic energy for small J/t. The exchange 
hole between parallel spins in the Gutzwiller wavefunction is consequence of the Fermi statistics. In the t — J model 
the motion of electrons with spin a is blocked also by the electrons with spin —a. A favourable kinetic energy is 
then obtained in a state in which an additional exchange hole in between electrons with opposite spins is simulated. 
This is exactly what the Jastrow factor defined by Eq. (|^) does, as it does contain short-ranged repulsion between 
all electrons, regardless of their relative spin orientation. Consequently the probability of finding two n.n. sites both 
occupied decreases with increasing 1/(2T) (for 5 = 1) and so does iV^o) as evidenced in Fig. At quarter filling 
there is one particle for every second site, which leads to a large increase of the probability of finding two n.n.n. 
sites both occupied and consequently Nt^^ increases with increasing 1/(2T). The extreme limit 1/(2T) — > oo is 
energetically not advantageous, as it corresponds to a charge-density wave with a delta function at (it, ir) and zero 
kinetic energy. Variationally we find the optimal value for 1/(2T) « 0.15 (within the constraint 5 = 1) for values of 
J < 0.25t. Within the accuracy of the variational calculations this agrees fairly well with the 1/(2T) w 0.2 which 
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would give the best result for Nq. 

The fact that the system wants to simulate an exchange hole between all electrons, regardless of the relative spin 
orientation has a natural explanation within the spinless Fermion model. In the case the holes would form a distinct 
Fermi-sea, as illustrated in Fig. (|l]), a exchange hole in between electrons of both spin orientations would be one of 
the consequences. Even a partial formation of a hole Fermi-sea, i.e. with a washed-out Fermi-surface, would also 
result in an effective, short-ranged charge-charge repulsion, as observed in the data for 7V q . Note that at quarter 
filling it/Hf ~ 1. We therefore interpret the data presented in Fig. (|^) as evidence for spin-charge separation at small 
lengthscales. 

In Fig. (^) we present the results at quarter filling for 7V q for the state with the best kinetic energy, allowing for 
any value of S. The optimal state is given by 1/(2T) « 0.15, S m 0.6, which is a Luttinger-liquid state J(|. It has 
qualitatively the same features as the states with n.n. correlations only, shown in Fig. (|J), but more washed out. 
In Fig. (|J) we present the results for A q for the Gutzwiller state for the cluster with L = 244 sites and Nh = 124 
and Nh = 52 holes, corresponding to electron densities n « 0.49 and « 0.79 respectively. The solid lines are the 
corresponding curves for spinless Fermions. We sec that near half- filling the magnitude of AT q , which again is very 
close to that of the spinless fermions ||, is reproduced by the Gutzwiller wavefunction. The data points for states 
1/(2T) > do not differ much from the data of the Gutzwiller state at these fillings, as density fluctuations are 
reduced by the projection operator. 

It is of interest to compare the variational Monte Carlo data for the Gutzwiller wavefunction with the predictions 
of the Gutzwiller approximation formula (GAF), which relates, by simple counting arguments p"4j ], the iV q in the 
projected state, as defined by Eq. (||), with the TVq ' in the unprojeced state, defined by 

JV<°> = l/ J L^e^( r *-^)(Vo|An i Ari J #o)/(Vo|V'o) (6) 

via the approximative formula 

N - - P> 

We have included in Fig. (Q) the predictions of the GAF, as given by Eq. (0) (dashed lines). The overall magnitude of 
iV q is well captured by the GAF, thought the variational Monte Carlo data has the tendency to be smoother at small 
dopings. The deviation of the data for the Gutzwiller state at n » 0.79 from the corresponding curve for spinless 
Fermions for small q-vectors is open to interpretation. Note, that no reliable results do exist for iV q in this region 
[^],^| . For the static spin-spin structure factor the Gutzwiller approximative fromula reads as 

In Fig. (||) we present Eq. (||) at quarter filling together with the variational Monte Carlo data for Sq for the Gutzwiller 
wavefunction (filled triangles) , the best Luttinger liquid wavefunction (filled circles) and the Fermi- liquid wavefunction 
with n.n. correlations only (filled squares). The overall magnitude is well reproduced by the GAF. Note that the 
Luttinger liquid state would have in the thermodynamic limit a cusp at 2kp with an algebraic singularity (at quarter 
filling 2kp = (-7T, 7r) along the (1,1) direction). 

In conclusion we have shown that the known features of the static charge-charge structure factor of the 2D t — J 
model are dominated by short-ranged correlations. We interpretat the results as evidence for spin-charge separation 
at small lengthscales. 
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FIG. 1. An illustration of the Brillouin zone of the 12 2 + 10 2 = 244 lattice with antiperiodic boundary conditions. The 
dashed lines indicate the location of the particle and the hole fermi surface at a filling n — (244 — 124)/244 ~ 0.49. The full 
circles/squares denote the k-states occupied by particles/holes on the finite lattice. Note that both the particle and the hole 
fermi-surfaces are closed-shell configurations. 

FIG. 2. The charge-charge structure factor, N q for various q-directions in the 2D Brillouin zone at quarter filling 
(n = 120/244). The filled circles are the results for the Fermi-liquid wavefunction with n.n. density-density correlation in 
the Jastrow prefactor of strenght 1/(2T) = 0.0, 0.15 0.2 and 0.25 respectively. The full line is the N q of spinless fermions, the 
dotted lines are guides to the eye. 

FIG. 3. The charge-charge structure factor, N q for various g-directions in the 2D Brillouin zone at quarter filling. The filled 
circles are the results of the optimal Luttinger-liquid wavefunction and the full line is the N q of spinless fermions, the dotted 
lines are guides to the eye. 

FIG. 4. The charge-charge structure factor, N q for various q-directions in the 2D Brillouin zone for fillings 
n = (244 - 124)/244 w 0.49 and n = (244 - 52)/244 « 0.79. The filled circles are the results for the Gutzwiller wave- 
function, the dashed line the predictions of the Gutzwiller approximation formula (GAF) and the full lines the respective N q 
of spinless fermions. The dotted lines are guides to the eye. 

FIG. 5. The spin-spin structure factor, S q , for various q-directions in the 2D Brillouin zone at quarter filling. The results 
for the Gutzwiller wavefunction (triangles), for the optimal Luttinger-liquid wavefunction (circles) and for the correlated 
Fermi-liquid wavefunctions (squares) are compared with the Gutzwiller approximative formula (GAF), given by the dashed 
line. The dotted lines are guides to the eye. 
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